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f— ^ has been previously considered intractable. 

o 

(N . 

q I PACS numbers: 02.50.Ey, 05.40Fb 

^! 1. Introduction 

The problem of a random walk on a two-dimensional lattice with a single interior source 
point and absorbing boundaries was first considered by Courant et al. |I| in 1928 where 
general properties of the solution were discussed. This problem was solved exactly in 
1940 § for the case of random walks on a square lattice with rectangular absorbing 
boundaries. The problem on a triangular lattice with finite absorbing boundaries 
was considered in 1963 where the exact solution was given for an approximation 
of the problem using straight boundaries rather than the true zig-zag boundaries of 
the triangular lattice. Indeed the authors of this paper remarked that "An explicit 
solution of the difference equation can hardly be obtained if these boundary conditions 
are used." Other variants of the problem on the square lattice have been solved exactly 
||. |4| however the problem on the triangular lattice with true zig-zag triangular lattice 
boundaries has remained unsolved. In this paper we give the exact solution for this 
problem. 

The problem of random walks on finite lattices with absorbing boundaries is 
fundamental to the theory of stochastic processes || and has numerous applications. 
These include potential theory 0, electrical networks pj, surface diffusion || and 
Diffusion-Limited Aggregation (DLA) ||. For example, the exact results for the square 
lattice problem have been used to expedite the growth of large DLA clusters on the 
square lattice ]l0| and the formulae derived in this paper could similarly be used on the 



triangular lattice. 



Random walks with absorbing boundaries 



2 



2. Field equations, boundary conditions and absorption probabilities 

A schematic illustration showing random walk pathways (dashed lines) on an equi- 
angular triangular lattice is shown in Fig. I. A natural co-ordinate system for the 
triangular lattice is to label the lattice vertices by the intersection points (p, q) of 
horizontal straight lines p — 0, 1, 2, . . . m+1 and slanted straight lines q — 0, 1, 2, . . . n+1 
- parallel to the p axis shown in Fig. 1. In this co-ordinate system the expectation that 
a random walk on the lattice visits a site (p, q), before exiting at a finite boundary, is 
given by the non-separable difference equation 

F(p, q) = i [F(p, q-l) + F(p, q + l) + F(p- 1, q) 

+F(p + l,q) + F(p + 1, q - 1) + F(p - 1, q + 1)] . 

In this paper we adopt a different co-ordinate system, which may at first appear less 
natural but it has the advantage that the governing equations are separable. 




q=0 1 2 3 4 5 6 7 q=8 

Figure 1. An equi-angular triangular lattice with natural (p, q) co-ordinates. 



In the co-ordinate system used here we label the vertices of the triangular lattice 
by the intersection points (p, q) of horizontal straight lines p = 0, 1, 2, . . . m + 1 and 
vertical zig-zag lines q = 0, 1,2, ...n+1 (see Fig. 2). In the absence of any source terms 
the expectation that a random walker visits a lattice site (p, q) is given by the coupled 
homogeneous difference equations 



F(p, q) = - [F(p - l,q - 1) + F(p,q - 1) + F(p + l,q - 1) 
+F(p +l,q) + F(p, q + l) + F(p-l,q) 

F(p, q) = l - [F(p -l,q) + F(p, q -l) + F(p+ 1, q) 

+F(p + 1, q + 1) + F(p, q+l) + F(p-l,q+l) 



(1) 



(2) 
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where F(p, q) denotes the field value at an even p co-ordinate and F(p, q) denotes the 
field value at an odd p co-ordinate. This separation of even and odd field equations, 
which is not necessary for the problem on the square lattice, is central to our solution 
below. 




Region I Region II 



Figure 2. Co-ordinate system, represented by dashed lines, adopted in this paper for 
random walks on the triangular lattice with true zig-zag lattice boundaries. In this 
example m — 7,n — 7. Nearest neighbour sites, highlighted with open circles, are 
shown surrounding a source point, highlighted by an open box, at a — 4, b = 4. The 
absorbing boundary sites are connected by solid lines. 



Now consider a source term at (a, b) and absorbing boundaries at p = 0, m + 1 and 
q = 0, n+1. In Section 3 we give the solution for the case where m > 1 is odd. The same 
methods can be used with slightly altered equations to obtain the solution for m even. 
The special case of the strip with m = 1 has been detailed elsewhere ||11|| . Following 
McCrea and Whipple [§] we construct separate solutions, Fi(p, q), Fj(p, q) forq<b and 
F n (p, q), Fji{p, q) for q > b. The absorbing boundary conditions are thus given by 



F{0,q) =0, 


(3) 


F(m+l,q) =0, 


(4) 


Fj(p,0) =0, 


(5) 


Fj(p,0) =0, 


(6) 


F n (p,n+1) = 0, 


(7) 


F H (p,n+l) = 0. 


(8) 



The ommision of a subscript / or II in the above indicates that the same equations are 
satisfied by both Fj and F n . 

At q = b we have the matching conditions 

Fj(p,b) = F n (p,b), (9) 
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F/(p,6) = F n (p,b). (10) 
Finally, with the inclusion of the source term, we have the coupled inhomogeneous field 
equations at q — b, 

6Fj(p, b) = 65 p , a ((a + l)mod2) 

+ [Fj(p -l,b-l) + F T (p, 6 - 1) + Fj(p + 1,6-1) 

+F J (p+l,6) + F // (p,6+l) + F / (p-l,6)] , (11) 
6Fj(p,b) = 65p i(1 (amod2) 

+ [Fj(p -1,6) + Fj(p, 6 - 1) + Fj{p + 1, 6) 

+F n (p +1,6+1) + F n (p, 6 + 1) + F n (p -1,6+1)]. (12) 

In Section III we present the solutions to Eqs. (1)— (12). 

The probabilities for the random walking particle to be absorbed at any specified 
point on one of the four boundaries, p — 0,p — m + l,q — 0,q — n + 1, are readily 
obtained by averaging over the nearest neighbour expectation values subject to the 
boundary conditions, Eqs. (3)-(8). For example, 

P{2k - 1, 0) = \ \Pr(2k -1,1) + Fj(2k -2,1) + F T (2k, 1) 
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*=1,..=£±. (13) 



3. Solution 



3.1. Homogeneous equations - Separation of variables 

We begin by solving the homogenous equations, Eqs. (1),(2), subject to the absorbing 
boundary conditions, Eqs. (3), (4). The homogenous field equations separate, on using 

F(p,q) = P(p)Q(q), (14) 

F(p,q) = P(p)Q(q), (15) 

into 

GQ(q)-Q(q+l)-Q(q-l) = F(p - 1) + P(p + 1) = 
Q(q-l) + Q(q) ' P(P) 

QQ(q) - Qjg + i) - Q(g - i) _ p(p - i) + Pip + i) n7] 

Q(q + 1) + Q(q) Pip) K ' 1 } 

where k and A are separation constants. The two coupled equations for P and P can 
be readily decoupled into separate equations 

P(2k + 2) + (2 - \K)P(2k) + P(2k - 2) = 0, (18) 

P(2k + 3) + (2 - \K)P(2k + 1) + P(2k - 1) = 0, (19) 

where k is an integer. This second order linear difference equation for P, (P) on a lattice 
of even (odd) integers has the solution 

P{2k) = Afi k + B^ k , 



Random walks with absorbing boundaries 5 
where 

-(2 - Xk) + J(2-Xk) 2 -A 

^= 2 ' 

and A and B are arbitrary. It follows that the solution of P(p) (with p even) 
corresponding to the boundary conditions, Eqs. (3), (4), is 

P(p) = csin(^^) (20) 

with 

Ak = 4cos 2 (^M. (21) 
The solution for P(p) (with p odd) now follows from Eq. (17) as 

P(p) = — cos(^— ) sm — - . (22) 
k m + 1 m + 1 

Equations (16) and (17) can also be decoupled for Q(q) and Q(q) resulting in the 
same fourth order linear difference equation in each case; 

Q(q + 4) - (12 + \K)Q(q + 3) + (38 - 2Xn)Q(q + 2) 

-(12 + \K)Q(q + 1) + Q(q) = 0. (23) 

The zeroes of the corresponding characteristic quartic polynomial are in reciprocal pairs 
leading to the solutions 

Q(q) = Cl e aq + c 2 e~ aq + c 3 e^ + c 4 eA (24) 
Q(q) = c ie aq + c 2 e~ aq + c 3 e Pq + (25) 



cosh a = 3 + ^ + ^\A Ak ) 2 + 32Ak ' ( 26 ) 
Xk 1 



where 



cosh/3 = 3 + — - -yj (Xk) 2 + 32Xk. (27) 

The coefficients ci, 02,03,04 may be chosen arbitrarily with the coefficients £1,62,63,04 
then determined by substituting Eqs. fl24D , ( |25"1) into Eq. (16) or Eq. (17). From Eq. 
(17) we obtain 

h = ci/(a), c 2 = c 2 f(-a), c 3 = c 3 /(/3), c 4 = c 4 f(-f3) (28) 

where 

/( ) + ^ (29) 

2(3 — coshu;) 

By combining Eqs. (14), (15), (20), (22), (24), (25), (28), (29) we can write the solutions 
to the field equations, Eqs. (|l|), @, that satisfy the boundary conditions, Eqs. (3), (4), 
in the form 



F(p,q;A,B,C,D) = J2 s ™(^TT 
JT{ \m + l 



x 



A je aq (3 - cosh a) + Bje _aff (3 - cosh a) 
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+Cje? q (3 - cosh/?) + D je -^(3 - cosh/5) 



6 

(30) 



F(p, q- A, B, C, D) = £ sin f^- 

Vm + 1 



cos 



7Ty 



m + 1 



x 



(31) 



A ie a,? (e a + 1) + Bje- a9 (e _a + 1) 
+CW + 1) + Dje-toie-P + 1)] , 

where A,-, 5.,, Cj, Dj are arbitrary. 

Before imposing the remaining absorbing boundary conditions, Eqs. (5)-(8), we 
represent the solutions in the two regions: / : q < b; II : q > b by 

Fi(p,q) = F(p,q;A,B,C,D), 
Fj(p,q) =F(p,q;A,B,C,D), 
F n (p, q)=F(p,q;A,B,C,D), 
F n (p, q) = F(p,q;A,B,C,D). 

The absorbing boundary conditions, Eqs. (5)-(8), can be used to eliminate C, D, C, D. 
After an appropriate re-normalization of the coefficients A, B, A, B this yields 



Fi( P ,q) = J2 sin 



j'=i 



_KJP_ 

m + 1 



x [A T (q, a, P)Aj + A T (q, -a, 0)Bj] , 



(32) 



where 



Flip, o) 



Fn(p, q) 



, njp \ f nj 
2 J sin I r I cos 



j'=i 
x 



m + 1/ Vm + 1 
[A I (q,a,p)A j + A I (q,-a,P)B J 



Fn(p, q) = £ sin f — ^r) 



x 



m + 1, 
a, + A 7/ (g, -a, (3)B 3 



™ . fnjp\ ( nj 

> sin cos _ 

^ Vm + iy Vm+1, 

x [i J7 (g, a, /?)!,■ + A n (q, -a, (3)B 3 



Aj(q,a,P) = (3 - cosh a) (3 - cosh/3)2 sinh/5e ag 
- (3 -cosh / 3)7(-a,- / 9)e /3,? 
+ (3-cosh/3)7(-a,/5)e- /3,? , 



(33) 



(34) 



(35) 



(36) 



i 7 (g, a, /?) = 2e Q,? (e a + 1)(3 - cosh (3) sinh /? 
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- 7 (-a, + + 7 (-a, /3)(e^ + l)e"^, (37) 



A n (q,a,0) = (3 - cosh a) (3 - cosh ,3)2 sinh/?e ag 

- e (Q - /3)(n+1) 7 (-a, -/?)(3 - cosh^e^ 

+ e (Q+/3)(n+1) 7 (-a, /?) (3 - cosh /3)e"^, (38) 



i 7/ (g,a,/3) = (e a + l)(3-cosh/3)2sinh/?e Q,? 

- e ( Q - /3 )( n+1 ) 7 (-a, -^(e' 3 + l)e^ 9 

+ e (Q+/3)(n+1) 7 (-a, /?) (e"" + l)eA (39) 

and the function 

j (a, (3) = 4 cosh a — 4 cosh /3 — (3 — cosh (3) sinh a 

- (3-cosha)sinh/5. (40) 

The coefficients defined in Eq. (36) which are shown as implicit functions of a and (3 are 
also functions of j via Eqs. (21), (26), (27). Equations (I32|)-(p5|), are general solutions to 
the coupled homogeneous field equations, Eqs. (1),(2), that satisfy all of the absorbing 
boundary conditions, Eqs. (3)-(8). 

3.2. Inhomogeneous equations - Matching conditions 

By using the two matching conditions, Eqs. (9), (10), the four arbitrary constants, 
A, B, A, B, can be reduced to two arbitrary constants, A, B say. The solutions in Region 
II can then be written as 

F H {p, q) = £ sin ( JE Tt) l A n(<l> ® A J + A n^ "«> P) B i\ > ( 41 ) 

Fn{p, q) = X/ sin ' r ' cos 



x 



m + U Vm + 1, 
[i^(g, a, 0)^- + i^(g, -a, /^-J , (42) 



where 



and 



A* n (q, a, (3) = A H (q, a, (3)T X + A n (q, -a, /3)T 2 , (43) 
A* n (q, a, (3) = A n (q, a, (3)T 1 + A n (q, -a, (3)T 2} (44) 



Fi = A n (b, -a, (3)A I {b, a, (3) - A n (b, -a, ff) Aj(6, a, (3) 
A n (b,a,(3)A n (b,-a,f3) - Aii(b,a,f3)A n (b,-a,(3y 

r2 = A H (b, a, (3)Aj(b, a, (3) - A n (b, a, (3)Aj(b, a, (3) 

A II (b,-a 1 f3)A n (b,a 1 (3)-A n (b,a,p)A II (b 1 -a,(3y 
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Finally the remaining arbitrary constants A, B are determined from the requirement 

that the solutions satisfy the coupled inhomogeneous field equations, Eqs. (11), (12). 
This step is facilitated using the identity 

Or, a — > sin sin . 47 

P ' m+lp[ \m + lj \m + lj V ; 

Explicitly we find: 
i) a even, 

4 = IK (48) 

B . = J?-^(J^-) m ■ (49) 

m + 1 Sm U + \) WmU) - Uj)T 2 (j) ' l4yj 



ii) a odd, 



and 



A > - m B '- (50) 

p _ 12 . / TTja \ T 3 (j) 

j ~ m + 1 sm U + 1) T 2 {j)Uj) - T^mu) ' 



TiO') = 2A 7 (6, -a, /3) + 2AJ 7 (6 + 1, -a, (3) + - 1, -a, /3) 

+ ij J (6 + l,-a,/?)-6A J (6,-a,/?), (52) 



T 2 {j) = 6Aj(b, a, (3) - 2Mb, a, (3) - 2A* n (b + 1, a, /3) 

-ij(6-l,a,/3) -i^(6+l,a,/3), (53) 

T 3 (j) = 6A 7 (j, 6) - 2 cos 2 (^-)i 7 (& - 1, a,/3) 

m + 1 

- 2 cos 2 (-^-) A 7 (6, a, (3) - Mb -l,a,P) 
m + 1 

-^(fe+l,a,/3), (54) 

T 4 (j) = 2 C os 2 (^-)i 7 (& - 1, -a, (3) 
m + 1 

+ 2 cos 2 (-^-)ij(&, -a, /?) + A 7 (6 - 1, -a, (3) 
m + 1 

+ AJ J (6 + l,-a,/?)-6A 7 (6,-a,/?). (55) 

Our solution for the triangular lattice site expectation values with absorbing 
boundary conditions and a source point is finally given by Eqs. (32), (33) in Region 
I and Eqs. (41), (42) in Region II. The relevant quantities appearing in these equations 
are defined through the series of equations, Eqs. (21),(26),(27),(36)-(40),(43)-(46),(48)- 
(50). 
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Table 1. Absorption probabilities at the boundaries on a triangular lattice with a 
source point enclosed by zig-zag lattice boundaries; (a) exact results using the formulae 
derived in this paper, (b) results using equations (1), (7a), (7b), (7c) derived in Keberle 
and Montet p] for the straight line boundary approximation. 



[A) 




(h) 
( >>) 




P(0,1) 


.012247 


P(2,0) 


.019609 


P(0,2) 


.035646 


P(4,0) 


.039949 


P(0,3) 


.054827 


^(6,0) 


.057904 


P(0,4) 


.063296 


P(8,0) 


.065834 


P(U, o) 


.UODDOD 


P(1U, UJ 




P(0,6) 


.039811 


P(12,0) 


.042817 


P(OJ) 


.020737 


P(14,0) 


.024583 


P(0,8) 


.005743 


P(16,0) 


.007903 


P(1,0) 


.026040 


P(l,l) 


.033624 


P(2,0) 


.013797 


P(0,2) 


.036214 


P(3,0) 


.069523 


P(l,3) 


.088151 


P(4,0) 


.021476 


P(0,4) 


.053888 


P(l,8) 


.005743 


P(17,l) 


.015276 


P(2,8) 


.040253 


P(16,2) 


.051612 


P(3,8) 


.015039 


P(17,3) 


.038562 


P(4,8) 


.059725 


P(16,4) 


.075927 



4. Example 

Consider the case of a triangular lattice with m = 7, n = 7 and a source at a = 4, b = 4 
(Fig. 2). In the approximation to this problem using straight edge boundaries and 
the (p, q) co-ordinate system of Keberle and Montet [[| this problem corresponds to; 
m = 16, n = 7, a = 8, b = 4. 

We have calculated expectation values at the nearest neighbour lattice sites around 
the source and absorption probabilities at the boundaries; (a) using our exact results 
above, and (b) using the results of Keberle and Montet for the approximation to 
the problem - their equations (7a), (7b), (7c). In this small model lattice system we 
found that the straight line boundary approximation provides reasonable results for 
expectation values near the source (accurate to within a few percent) but provides poor 
results for the absorption probabilities (Table 1). Note that the 'absorption probabilities' 
for the zig-zag boundary co-ordinates which are inside the straight line boundary of 
Keberle and Montet || are not true probabilities in their solution. 

5. Conclusion 

Although our solution is somewhat unwieldly, we have nevertheless solved the underlying 
field equations for random walks on a finite triangular lattice with a single interior 
source point and zig-zag absorbing boundaries and have thus calculated the associated 
absorption probabilities. This problem was previously considered to be intractable 0. 
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We hope that our result will inspire further work in this area. 
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